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Abstract 

In this paper, we introduce a new class of bivariate distributions by compounding the 
bivariate generalized exponential and power-series distributions. This new class contains 
some new sub-models such as the bivariate generalized exponential distribution, the bivari¬ 
ate generalized exponential-poisson, -logarithmic, -binomial and -negative binomial distri¬ 
butions. We derive different properties of the new class of distributions. The EM algorithm 
is used to determine the maximum likelihood estimates of the parameters. We illustrate 
the usefulness of the new distributions by means of an application to a real data set. 
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1 Introduction 


The modeling of lifetime is an important aspect of statistical work in a variety of scientihc 
and technological helds. In this area, much new univariate distributions have been studied in 
the statistical literature recently. Interestingly, not much work has been done on the bivariate 
distributions mainly due to its analytical intractability. 

The two-parameter generalized exponential (GE) distribution has been introduced by 


Gupta and Kundul (jl999l i and it has the following cumulative distribution function (cdf) and 


probability density function (pdf), respectively: 

FGEix]a,X) = {1 — e~^^)°‘, x > 0, a,A>0, 

fGE{x;a, X) = aXe~^^{l — e~^^)°‘ x > 0, a, A>0. 


( 1 ) 

(2) 
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The hazard function of GE distri bution can be increasing, dec reasing and constant, but it 


cannot be bathtub shaped. Therefore, 


Mahmoudi and Jafaril (j2012l i introduced the generalized 


exponential-power series (GETS) distributions by compounding GE distribution with power se¬ 


ries c lass of distributions. The proposed class includes GE, Poisson -exponential (|Gancho et al. 


power series ( 
Recently, 


Roresetal 


201lh . compl ementary exponenti al-geometric (jLouzada et al.l . l201ll ). and complementary exponential- 


201311 distrib utions. 

Kundu and Guota ( 2 OO 9 I 1 extended the GE distribution to a bivariate distribu¬ 


tion. The cdf of the bivariate generalized exponential (BGE) with parameters ai, 0 : 2 , 03 and 
A is given by 


Fbge {xi,X2 ‘,«!, a2, as, A) = 


(1 — e _ g ^X2'^a2 Xl < X2 


(1 


_ _ f,-\x2\a2+a3 


if Xi > X2- 


( 3 ) 


Note that the BGE distribution has both an absolutely continuous part and a sing ular part 


similar to the bivariate expon ential distribution reported in 


bivariate models proposed by 


Sarhan and BalakrishnanI (120071 ) 


Marshall and Olkii] (j 19671 ) and the 


In this paper, we compound the BGE distribution and power series class of distributions 
and define a new class of bivariate distributions. It contains the BGE and GEPS distributions 
and is called the bivariate generalized exponential-power series (BGEPS) distributions. This 
paper is organized as follows. In section [2l we introduce the BGEPS model and obtain some 
properties of this new family. Some special models are studied in detail in Section[3l We propose 
an EM algorithm to estimate the model parameters in Section |4l A real data application of 
the BGEPS distributions is illustrated in Section [5l 


2 The BGEPS class 

A random variable N follows the power series distribution if it has the following probability 
mass function 


P{N = n) = 


c{ey 


n = 1,2, 


( 4 ) 


where > 0 depends only on n, C{9) = 9 G (0, s) (s can be 00 ) is such 

that C{9) is finite. Table [U lists some particular cases of the truncated (at zero) power se¬ 
ries distributions (geometric, Poisson, logarithmic, bin omial a n d neg ative binomial). Detailed 


Noack (195C). Here, C'{9), C"{9) and 


properties of power series distribution can be found in 
C'"{9) denote the first, second and third derivatives of C{9) with respect to 9, respectively. 

Now, suppose {(Xin,X 2 n)','n = 1,2... } is a sequence of independent and identically dis¬ 
tributed (i.i.d.) non-negative bivariate random variables with common joint distribution func- 
































Table 1: Useful quantities for some power series distributions. 


Distribution 

an 

C{9) 


C'{9) 

C"{9) 

C'"{9) 

S 

Geometric 

1 


-1 

(1-0)-^ 


6(1-0)-^ 

1 

Poisson 

n\~^ 

e^-l 



e® 


00 

Logarithmic 

n~^ 

- log(l - 

9) 

(l-0)-i 

(l-0)-2 

2(1 - 0)-3 

1 

Binomial 

0 

(i + oy- 

- 1 

k 

fc(fc-i) 

k{k-l){k-2) 

(6»+l)3-'= 

00 

Negative Binomial 

/n-l\ 

gk 



fc(fc-|-26»-l) 


1 

U-J 



(l_0)fe+i 

g2-k^l_g)k+2 

03-fc(l_0)fc+3 



tion Fxi-, •), where X = {Xi,X 2 )'■ Take to be a power series random variable independent 
of {Xu,X 2 i). Let 

Tj = maxjXii,... XjAr}, i = l,2. 


Therefore, for the joint random variables Y = (Yi, 12 )^ we have 


P{Yi < 2/1,12 < y2, N = n) = FYj^x2\N{yi,y2\n)P{N = n) = (Tx(2/i, 2/2)) 
Therefore, the joint cdf of = (Ti, 12 )^ becomes 

FY{yuy.} = ■ 

In this case, we call Y has a bivariate F-power series (BFPS) distribution. 

The corresponding marginal distribution function of Ij is 

C{0FxM)) 


c{ey 


Pviiyi) = 


m 


i = 1,2. 


(5) 


( 6 ) 


In recent years many anthors have conside red this univariate class: for example the GEPS 


distribntion by 
distribntion by 


Mahmondi and .lafaril ([2013) and the complementary exponential-power series 


Flores et al 


(1201^ ) among others. 


Remark 2.1. If we consider Zi = min {Xji,... i = 1,2, another class of bivariate 

distribution is obtained with the following joint cumulative survival function: 


FzYZ 2 {yi,y 2 ) = P{Zi > yi,Z 2 > 2 / 2 ) = 


C{eFxiyi,y2)) 

0(9) 


where Fx{yi,y 2 ) = P{Xi > yi,X 2 > 2 / 2 )- Several papers have stu died the univariate case 
of this class: for example the exponential-power seri es distribution bvlChahkandi and Ganial 


(200m) and the Weibutll-power series distribution bv lMorais and Barreto- Souza l20li) among 
others. 










































In this paper we take F to be the bivariate generalized exponential given in ([T]). Therefore, 
we consider the bivariate generalized exponential-power series (BGEPS) class of distributions 
which is defined by the following cdf: 


FY{yi,y2) = < 


cm 

C(0) 

C {OFce {yi ;«! +«3 , (j /2 ;«2 , a)) 

c{e) 

C{SFaE{yi\ai,>~)FGE(y2\a2+0‘3,>')) 

c{e) 


if yi < 2/2 

if 2/1 > 2 / 2 , 

if 2/1 < 2/2 

if 2/1 > 2 / 2 , 


( 7 ) 


We denote it by BGEPS(ai, a 2 , as, A, 0). 


Proposition 1. Let FY{yi,y2) be the cdf of BGEPS distributions given in d?]). Then 

OO 

TV (2/1,2/2) = '^PnFBGE{yi,y2;nai,na2,na3,X), 

n=l 

where pn = P{N = n) = . 

Proposition 2 . Let (Ei, 12)/o/^ows BGEPS (ai, 02,03, A, 0 ) distribution. Then 

1 . Eaeh Yi has a GEPS distribution with parameters a* + 03, A and 9 . 

2 . The random variable U = max(Ei,y2) has a GEPS distribution with parameters 01 + 02 + 03, 
A and 9 . 


3 . If C ( 9 ) = 9 , then Y has a BGE distribution with parameters oi, 02, 03 and A. 


4 . P(Ei <E 2 ) = 


ai 


01-1-02-1-03 ' 


Theorem 2 . 1 . Let Y /las a BGEPS (oi, 02, 03, A, d) distribution. Then the joint pdf of Y is 

{ /i (2/1,2/2) if 0<yi<y2 

/ 2 ( 2 /i, 2 / 2 ) if 0 < 1/2 < 2/1 ( 8 ) 

foiy) if 0 <yi = y 2 = y, 

where 

Q 

/i (2/1,2/2) = ^7 ^/ge(2/i;«i+ «3,A)/GE(2/2;a2,A) [6'FG£;(2/i;ai+«3,A) 

X-Pge( 2 / 2 ; 0 L 2 , X)C'' { 9 FGEiyi] ai + 03, X)FGEiy 2 ', 0 ^ 2 , A)) 


/2 (2/1, 2/2) 


+C' { 9 FGE{yi] ai + 03 , A)Egs( 2 / 2 ; 02, A))], 

Q 

qJ^J'ce (2/1; «i. A) fGE (2/2; a2 + as. A) [ 0 Fge( 2 /i; ai, A) 
xFGE{y 2 ', OL2 + 03, A)C" ( 0 Fg£;( 2 /i; ai, A)Fge(2/2; 02 + 03, A)) 


(9) 


/o (2/) = 


+C' {9FGE{yT, ai, X)FGE{y2', a2 + 03 , A))], 
9a3 


( 10 ) 


C ( 9 ) (oi + 02 + 03 ) 


fGE (2/; oi + 02 + 03 , A) C' { 9 Fge (2/; oi + 02 + 03 , A)). 


( 11 ) 












Proof. It is obvious. 


□ 


As a special case, consider C{6) = 9 + 6^^. It is also considered by 


Mahmoudi and Jafari 


(j2012l l. For A = 1 and other values of the parameters, the pdf of the BGEPS class of distribu¬ 


tions are depicted in Figured] 



ti.14- 





Figure 1: The pdf of the BGEPS class of distribution for some values of parameters: ai = 
02 = as = 1, 0 = 1 (left top), ai = 02 = as = 1, 0 = 2 (right top), oi = 02 = as = 1, 0 = 0.5 
(left bottom), oi = 02 = as = 10 , 0 = 1 (right bottom). 


Remark 2.2. Since 9C" (0) + C (0) = ^ o,nd C (0) = ^ , we have 

00 

/i (01,2/2) = '^PnfcE {ypnai +na 3 ,X)fGE {y2]na2,X), 

n=l 

00 

72 ( 01 , 02 ) = '^PnfcE {yr,nai,X) fcE ( 02 ;nQ !2 + nQ!s,A), 


/o (0) = 


n=l 

as 


ai + a2 + as 


'^PnfcE iy,nai + na2 + na3,X), 


n=l 


where pn = P {N = n) = and fcE (■; 11 a, A) is the pdf of GE distribution with parameters 










na and A. Note that /ge(-; na, A) is the pdf of random variable max(?7i,..., Un) where Ui’s 
are independent random variables from a GE distribution with parameters a and A. 


Corollary 2.1. The joint pdf of the BGEPS distributions provided in Theorem \2.1\ can be 


written as 


where 


fY{yi,y2) = 


Oi + 0i2 

Oi + 02 + 0:3 


ga{yi,y2) + 


03 


Ol + 02 + «3 


gs{y), 


( 12 ) 


1 , 2 / 2 ) 


Ol + O2 + 03 
Oi + O 2 


/i( 2 / 1 , 2 / 2 ) if 2/1 < 2/2 

/ 2 ( 2 /i, 2 / 2 ) if 2/2 < 2 / 1 , 


gs{y) 


e 

W) 


fcE ( 2 /; Ol + 02 + 03, a) C' (OFgE ( 2 /; Ol + 02 + 03, A)) 


if 


2/1 = 2/2 = y, 


and 0 otherwise. Glearly, gai-,-) is the absolute continuous part and gs{.) is the singular part. 
U ct 3 = 0 , does not have any singular part and it becomes an absolute continuous density 
function. Note that gs{.) is the pdf of GEPS distribution with parameters oi + 02 + 03 and A. 


Proposition 3. The conditional distribution of Yf given Y 2 < 2/2 is an absolute continuous 


distribution function with the following cdf: 


P{Yi<yi\Y2<y2) 


)“i +“3 (i_e-^ 2 ' 2 )“ 2 ) 
C(6»(l-e-^a2)“2+“3) 
C(e(l-e -^!'2 )“ 2+“3 (l-e-'i^'l )“l) 
C( 6 »(l-e-^« 2 )“ 2 +“ 3 ) 


if 2/1 < 2/2 
if 2/2 < 2 / 1 - 


Proposition 4. The limiting distribution of BGEPS when 0 ^ O'*’ is 


lim iV ( 2 / 1 , 2/2) 
e-)-o+ 


0 ^ 0 + C {9) 0 ^ 0 + Z)n=l 

ac{Fx {yi,y 2 )r + EZc+i^n9^-'^{Fx ( 2 / 1 , 2/2))" 

0 ™+ «c + E“=c+l«n0"-" 

{Fx {yi,y2)T 

I (1 _ (1 _ ^-Ay,y(»2+as) ^ 


which is the pdf of a BGE distribution with parameters cai, 002 , C 03 and X, where c = min{n G 
N : a„ > 0}. 


For the joint random variables {Yi,Y 2 , N), consider equation ([5|) when X has a BGE 
distribution. Since (Yi,y 2 |-iY = n) has a BGE with parameters noi, no 2 , 1103 , and A, the joint 
pdf of iYi,Y 2 ,N) is 


fYi,Y2,N {yi,y2,n) 


^^finiyi,y 2 ) if 2/1 < 2/2 

< ^/2n(2/l,2/2) if 2/2 < 2/1 

v^foniy) if 2/1 = 2/2 = y, 















where 


fin{yi,y2) = n^X\a, + 

f2n{yi,y2) = n^\\cx2 + as)aie-^y^-^y^{l - 

foniy) = nXa3e-^y{l - ^-^y)n(a,+a,+as)-l _ 


The conditional probability mass function of N given Yf = yi and I2 = ?/2 is 


fN\Yi,Y2 (^11/1)1/2) 


< 


n'^an{dKi{yi,y2))"~^ 

hi{yi,y2) 

n^an{SK2{yi,y2))"~^ 

h2{yi,y2) 

na„{dKo{y))"~^ 

hoiy) 


if 2/1 < 2/2 

if 2/2 < 2/1 

if 2/1 = 2/2 = y, 


where 


^1(2/!) 2/2) = OFge (2/1; ai + A) Fge (1/2; a2, A) C" {OFge (2/1; «i + «3, A) Fg_b (2/2; 02, A)) 
+C' {OFge (2/1; «i + 03, A) Fg_b (1/2; a2, A)), 

^1(2/!) 2/2) = (yi; ai, A) Fg_b (^2; a2 + 0:3, A) C" {OFge (2/1; «1) A) Fge (2/2; a2 + a3, A)) 

+C' {OFge (2/1; «1) A) Fg_b (^2; a2 + a3) A)), 

^ 0 ( 2 /) = C''(6'Fg£;(2/;«i+ a2 + a3,A)), 

and 


Ki (2/1,2/2) = (1 - _ e-Aj;2)«2 = (2/2; «2, A), 

K2 (2/1,2/2) = (1 - - e-^^2)«2+a3 = (yi; ai, A) Fgs (^2; 02 + «3, A), 

Ko{y) = {l-e-^yr^+^^+'^^=FGE{y;ai + a2 + a3,X). 


Since O^C” { 0 ) + 300 ” { 0 ) + C { 0 ) = n^an 0 ^-\ OC { 0 ) + C { 0 ) = n'^UnO^-^ and 

F' { 0 ) = nan 0 ^~^-, therefore, we can obtain the conditional expectation of N given Yi = yi 

and y2 = 2/2 as 


E{N\yi,y2) 


< 


Ri{yi,y2) 

hi(yi,y 2 ) 

R2{yi,y2) 

h2{yi,y2) 

eKo(y)C''(eKo(y))+C'(eKo(y)) 

ho(y) 


if 2/1 < 2/2 

if 2/2 > 2/1 

if 2/1 = 2/2 = y, 


(13) 


where 


Ri{yi,y2) = {0Ki{yi,y2)fC'"{0Ki{yi,y2)) + 36»Fri(yi, y2)C'"(6'Fi(yi, ^ 2 )) 


+C'{0Ki{yGy2)), i = l,2. 



3 Special Cases 


In this section, we consider some special cases of BGEPS distributions. 


3.1 Bivariate generalized exponential-geometric distribution 


When C (6) = {0 < 9 < 1), the power series distribution becomes the geometric distribution 

(truncated at zero). Therefore, the cdf of bivariate generalized exponential-geometric (BGEG) 
distribution is given by 


Fy ( 2 / 1 , 2 / 2 ) 


(l_0)(l_e-A!/l)“l(l_e-A!/2)“2+“3 

l_0(l_e->w)“i(l-e-^!/2)“2+“3 


if 2/1 < 2/2 
if 2/1 > 2/2, 


and its pdf is given in ([8|) with 


/i ( 2 / 1 , 2 / 2 ) 
72 ( 2 / 1 , 2 / 2 ) 


(1 - ^)/g£;(2/i; «i + «3, A)/ge(2/2; 02, A) 
(1 - ^)/ge(2/i; «1, A)/ge(2/2; 02 + 03, A) 


1 + 6 IFge( 2 /i; oi + 03 , A)Fge( 2 / 2 ; 02 , A) 

(1 - 9FGE{yi] Oi + 03, X)FGE{y2',CX2, A ))3 

1 + 6 IFge( 2 /i; oi, A)Fge( 2 / 2 ; 02 + 03 , A) 

(1 - OFcEiyi] oi, A)Fge(2/2; 02 + 03, A))^ 


fo{y) 


_ (1 - 9)a3fGE{y; oi + 02 + 03 , A)_ 

(ai + a2 + 03)(1 - 9 FGE{y,ai + 02 + 03, A))^ 


Remark 3.1. When 9* = 1 — 9, we have 


Fy ( 2 / 1 , 2 / 2 ) = < 

It is also a cdf for all 9* > 0 (see 




l-(l-6»*)(l-e-^2'l)“l 3(']^_g-A!,2)“2 

0*(l_e-A!,i)“ljl_g-A«2)“2+«3 

l- (l- 6 >*)(l-e-^^'l)°^(l-e - ^ 2)°'^ ''~°^3 


Marshall and Olkin. 


1997). In 


if 2/1 < 2/2 
if 2/1 > 2/2- 

fact, this is in Marshal-Olkin 


bivariate class of distributions. 


3.2 Bivariate generalized exponential-Poisson distribution 

When o,n = ^ and C (9) = e^ — 1 {9 > 0), the power series distribution becomes the Poisson 
distribution (truncated at zero). Therefore, the cdf of bivariate generalized exponential- Poisson 
(BGEP) distribution is given by 


Fy ( 2 / 1 , 2 / 2 ) 


)“l+“3 (•i_,-Ab2 )“2 
e^'-l 

e^-l 


and its pdf is given in ([8|) with 


if 2/1 < 2/2 
if 2/1 > 2/2, 


71 ( 2 / 1 , 2 / 2 ) = 6'7ge(2/i;oi- h03,A)7GE(2/2;02,A)e' 


fiFoE (yi ;ai -I-03 MFge (y2 ;«2, A) — 1 

















/2 (yi,y2) 


/o {y) 


X [0FGE{yi]ai + 03 , X)FGE{y2] 02, A) + 1] , 

OfcE (yi; ai, A) Jge (y2; a2 + as, A) e®''^GB(yi;ai,A)FGB(j/2;a2+a3,A)-i 

X [OFge (yi; ai, A) Fge {y2\ot2 + as, A) + l], 

7 -- t/gf (y; ai + 02 + as, A) e^FGE{y,»i+c. 2 +»z,X)-i_ 

(ai + 02 + asj 


3.3 Bivariate generalized exponential-binomial distribntion 

When an = and C (6) = {9 + 1)^ — 1 {6 > 0), where k{n < k) is the number of replicas, 
the power series distribution becomes the binomial distribution (truncated at zero). Therefore, 
the cdf of bivariate generalized exponential- binomial (BGEB) distribution is given by 


Fy (yi,y2) = < 




IT-; - if yi < y 2 


( 0 + 1 )'=-! 


if yi > y 2 . 


and its pdf is given in 

/i(yi,y 2 ) = 


with 

ke 


/2(yi,y2) — 


u -/GE(yi;ai-Fa3,A)/GF(y2;a2,A) 

(6» + if -1 

X [9Fge (yi; ai -l- as. A) Fge {y 2 ‘, a2, A) -|-1]^ ^ 

X [k6FGE{yi',ai + as, \)FGE{y2]oi2, A) -|- 1] , 

kO 

-fcE (yi; ai. A) /ge (y 2 ; a 2 -I- as. A) 


/o(y) = 


{e + if - 1 ' 

X [OFge (yi; ai. A) Fge (y2; a2 -I- as. A) -|- 1 ]^ ^ 

X [k6FGE{yi]Oii,\)FGE{y2',Oi2+a3A) + 1] , 

kOasfGEiy, ai+a2 + as. A) 


{9FGE{y, ai + a2 + as, A) 1) 


k-l 


((0 -|- 1) “ l)(ai -|- a2 -|- as) 

3.4 Bivariate generalized exponential-logarithmic distribution 


When On = A and C (0) = —log(l — 6) {0 < 9 < 1), the power series distribution becomes 
the logarithmic distribution (truncated at zero). Therefore, the cdf of bivariate generalized 
exponential- logarithmic (BGEL) distribution is given by 


•Fv(yi,y2) = 


log(l-0) 

log(l-0(l-e-^2'l)“l(l-e“^!'2)“2+“3) 

log(l-0) 


if yi < y 2 
if yi > y 2 , 


and its pdf is given in ([8]) with 


_ -OfoE (yi; ai + as, A) Jge {yy, a 2 , A)_ 

log(l — 9){l — OFge (yi; ai -|- as. A) Fge (y 2 ; a 2 , A))^ 


fi (yi,y2) 












/2 ( 2 / 1 , 2 / 2 ) 


/o (y) 


_ -(^fcE fcE (2/2; «2 + « 3 , A)_ 

log(l — 0)(1 — 6 Fce (2/1; «1) ■^) -Pge (2/2; 02 + as, A))^ 
_ —Oa^fcE jy, «! + «2 + Q 3 ; A)_ 

log(l - 6 ») (ai + 02 + as) (1 - OFge (y, ai + 02 + as, A))' 


3.5 Bivariate generalized exponential- negative binomial distribution 


When an = (fc_i) and C{9) = (j^)^ {0 < 9 < 1), the power series distribution becomes the 
negative binomial distribution (truncated at zero). Therefore, the cdf of bivariate generalized 


exponential- negative binomial (BGENB) distribution is given by 


( 2 / 1 , 2 / 2 ) 


(l_0(l_e-Aw)“l+“3(i_g-Aj,2)“2)'' 
, (l-e(l-e-^W)“i(l-e-^!' 2 )“ 2 +- 3 )'' 


and its pdf is given in ([8|) with 


if 2/1 < 2/2 
if 2/1 > 2/2, 


/i ( 2 / 1 , 2 / 2 ) 


72 ( 2 / 1 , 2 / 2 ) 


fo{y) 


k{l - 9)’^fcEivi^ai + 03 , A)/Gg(?/2; 02, A) 
(1 — 0 Fge(2/i; ai -I- as, X)FGE{y2] 0 : 2 , A))^"*"^ 
X [A: -I- 0Fg_b(2/i; ai + as, A)Fg'e( 2/2; a2. A)] , 


71 /c —1 

G£; 


( 2 / 1 ; ai -I- as, A)Fg^^ 


( 2 / 2 ; a 2 . A) 


A:(l - 6')^/GE(2/i;ai, A)/GE(2/2;a2-h as. A) fc_i, \\z7fc-i 

--T—--77---7- T^^ge (2/i;ai,A)Fcs 

(1 — 9FgeKVi] Oil, \)FGE\y2', a 2 + as, A)) 

X [fc -I- 9FGE{yi]Oii, X)FGE{y2', 02 + as. A)], 

fcas(l - 9)^fGE{y]Oii + 02 + 03, A)F^7n2/;«i + «2 + 03 , A) 

(oi -|- 02 + a 3 )(l — 9Fge(jJ'', Oil -|- 02 -I- 03 , A)) 


( 2 / 2 ; 02 + 03 , A) 


4 Estimation 


In this section, we consider the estimation of the unknown parameters of the BGEPS dis¬ 
tributions. Let ( 2 / 11 , 2 / 12 ) ,•••, ( 2 /mi, 2 /m 2 ) be an observed sample with size m from BGEPS 
distributions with parameters 0 = (oi, 02 , 03 , A, 0)^ Also, consider 

h = {i- yii = 2/2* = 2/i} , h = {i- 2/1* < 2/2*} , ^2 = 1 + yu > 1/2*} , 


and 

mo = |/o|, mi = \Ii\, m2 = 1/2!, m = mo + mi + m2. 

Therefore, the log-likelihood function can be written as 

log (/o (2/*)) + ^ log (/i (2/1*, 2/2*)) + X/ (2/1*, 2/2*)), (14) 

2€/o i £ l 2 

where fi, /2 and /o are given in Q, (fTO|) and respectively. We can obtain the MLE’s 
of the parameters by maximizing i ( 0 ) in (I14|) with respect to the unknown parameters. This 










is clearly a five-dimensional problem. However, no explicit expressions are available for the 
MLE’s. We need to solve five non-linear equations simultaneously, which may not be very 
simple. Th e maximization can be performed using a command like the nlminh routine in the 


R software (jR Development Core Teaml . 


2014 1. But, it is related to initial guesses. Theref ore, 


we present an expectation-maximization (EM) algorithm similar to lKundu and Devi (j2009l l to 
find the MLE’s of parameters. 

For given n, consider that independent random variables {Zi\N = n}, i = 1,2,3 have the 
GE distribution with parameters na* and A. It is well-known that 


{YilAT = n} = {max(Zi, Z 3 )|iV = n}, 


{Y 2 \N = n} = {max(Z 2 , Z 3 )\N = re}. 


Assumed that for the bivariate random vector (Yi,!^), there is an associated random 
vectors 

Ai = I ? and A 2 = 


0 ¥2 = Zi 

1 Y2 = Z^. 


Note that if Yi = I 2 , then Ai = A 2 = 0. But if 1} < Y 2 or Yi > Y 2 , then (Ai, A 2 ) is missing. 
If (Yi, Y 2 ) S Ii then the possible values of (Ai,A 2 ) are (1,0) or (1,1), and If (Yi,Y 2 ) E I 2 then 
the possible values of (Ai, A 2 ) are (0,1) or (1,1) with non-zero probabilities. 

We form the conditional ‘pseudo’ log-likelihood function, conditioning on N^ and then re¬ 
place N by E(A^|Yi, I 2 ). In the E-step of the EM-algorithm, we treat it as complete observation 
when they belong to Iq. If the observation belong to /i, we form the ‘pseudo’ log-likelihood 
function by fractioning (^ 1 ,^ 2 ) to two partially complete ‘pseudo’ observations of the form 
(yi,y 2 ,rei (0)) and {yi,y 2 ,U 2 (0)), where rei (0) and U 2 (0) are the conditional probabilities 
that (Ai,A 2 ) takes values (1,0) and (1,1), respectively. Since 


P(Z3 < Zi < Z2\N = n) = 
P{Zi < Z 3 < Z2\N = n) = 


Oi\(y.2 


(ofi -|- 03) [ai -|- 0:2 -|- 0:3) 
a2«3 

(cti -|- 03) [ai -|- Q!2 -|- U3) 


therefore. 


ui ( 0 ) = 


ai 


U2 ( 0 ) = 


as 


(15) 


ai -|- as Oil Y Q^3 

Similarly, If the observation belong to I 2 , we form the ‘pseudo’ log-likelihood function 
of the from {yi,y 2 ,vi (0)) and (yi, 2 / 2 , 1^2 (0)); where vi (0) and V 2 (0) are the conditional 
probabilities that (Ai,A 2 ) takes values (0,1) and (1,1), respectively. Therefore, 


a2 


a3 


a2 -I- a3' 


02 + «3 


Vi ( 0 ) 


V2 ( 0 ) 


(16) 














For brevity, we write ui (0), U 2 (©), vi (0), V 2 (0) as rti, U 2 , vi, V 2 , respectively. 

E-step: Consider bi = E{N\yii,y 2 i,&). The log-likelihood function without the additive 
constant can be written as follows: 

m 

4seudo(®) = log( 6 >) y^^bj- mlog{C{6)) + {uiq + 2mi + 2 m 2 )log(A) 

i=l 

+ {miui + m 2 ) log(Q;i) -I- (mi -|- m 2 Ui)log(a 2 ) 


+ {mQ + miU2 + m2V2)\og{a:i) - \C^yi+ ^ {yii + y2i)) (17) 

iG/o iG/iU /2 

+aiC^biQ{yi) + ^ &i<3(yii)) + a2(^ ftiQ(?/i) + ^ biQ{y2i)) 

jg/o iehuh i&Io iG/iU/2 

+a:^C^biQ{yi) + ^ biQ{yii) + ^ biQ{y2i)) 

2G/o 2€/i i£l2 

ig/o jg/lU/2 ig/iU/2 

where Q{y) = log(l — e~^y). 

M-step: At this step, ■^pseudo(0) is maximized with respect to ai, 0 : 2 , 03 , A and 9. For fixed 
A, the maximization occurs at 


di(A) = 

«2(A) = 

«3(A) = 


-{miui + m 2 ) 


Sig/o biQiyi) + Y2i£liUl2 biQiyii) 

— (mi -I- m 2 Vi) 
biQiyi) + X)iG7iU/2 biQ{y 2 i) 
-(mo + miU 2 -I- m 2 V 2 ) 


Sig/o biQiyi) + Ylieli biQiyii) + X)iG /2 biQiy 2 i 

and solving the following non-linear equation with respect to 9: 


(18) 

(19) 

( 20 ) 


Remark 4.1. When C (0) 


9C'i9) YZih 
Ci9) m 


( 21 ) 


then the solution of equation in (12X1) is 9 = 1 — 




Remark 4.2. IFe do not need to solve the equation in m, when C i9) = 9. In fact, the 
BGEPS distribution reduces to the BGE distribution. 


Finally, A can be obtained as a solution of the following equation: 

mo -|- 2mi -|- 2m2 


g(A) 


= A, 


(22) 


where 


g(A) = '^yi+ iyii + y2i) - ^ ikai -h 6 * 0:2 6 * 0:3 - 1 ) 


iGlo 2G/iU/2 


ieio 


Vie 

1 — e~^y* 










- 2 ^ ihai + has - 1 ) 

i&h 


Y1 (^*"2 - 1 ) 

ieh 


y2ie 

1 — e~^y^i 


- {hai - 1) 


yue 

1 — Q~^yii 


Y (^i«2 + has - 1) 


y2ie 

1 — ^~^y2i 


i£l2 iGl2 

The following steps can be used to compute the MLE’s of the parameters via the EM 
algorithm: 

Step 1: Take some initial value of 0, say 0*^°^ = q: 2 °\ 03 °^ )'. 

Step 2: compute h = E{N\yii,y 2 i]&^°^) 

Step 3: Compute ui, U 2 , vi, and V 2 - 

Step 4: Find A by solving the equation (l2^ . say A^^^. 

Step 5: Compute = Q!i(A(^^), i = 1,2,3 from (fT 8 ]l - (l 20 ]l . 

Step 6 : Find 9 by solving the equation (1^ . say 9^^\ 

Step 7: Replace 0^*^^ by 0^^^ = A’^^^ 0^^^), go back to step 1 and continue the 

process until convergence take place. 


,000000000 


n ^ I I I r 

5 10 15 20 25 30 



k 
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Figure 2: The log-likelihood values of BGEB (left) and BGENB (right) distributions for dif¬ 
ferent values of k. 


5 A real example 


MeintanisI (j2007l ) and is obtained from th e group stage of the UEFA 


The data set is given from 

Cha mpion’s League for the y ears 2004-05 and 2005-2006. In addition. 


and 


Kundu and Guotal (2009|) 


Kundu and Devi (j2009l l analyzed this data by using BCE and Marshall-Olkin bivariate 


Weibull distributions, respectively. The data represent the football (soccer) data where at 

















least one goal scored by the home team and at least one goal scored directly from a kick goal 
(like penalty kick, foul kick or any other direct kick) by any team have been considered. Here 
Yi represents the time in minutes of the first kick goal scored by any team and I 2 represents 
the first goal of any type scored by the home team. 

We divided all the data by 100. Then six special cases of BGEPS distributions are consid¬ 
ered; BGE, BGEG, BGEP, BGEB, BGENB, and BGEL. Using the proposed EM algorithm, 
these models are fitted to the bivariate data set, and the MLE’s and their corresponding log- 
likelihood values are calculated. The standard errors (se) based on the observed information 
matrix are obtained. The results are given in Table [2j 


Table 2: The MLE’s, log-likelihood, AIG, AICC, BIG, K-S, and LRT statistics for six sub¬ 
models of BGEPS distribution. 



Distribution 

Statistic 

BGE 

BGEG 

BGEP 

BGEB(A: = 30) 

BGENB(A: = 2) 

BGEL 

Ctl 

1.4452 

0.9964 

0.5644 

0.5980 

0.2538 

1.1871 

(s.e.) 

(0.4160) 

(0.4938) 

(0.5758) 

(0.5442) 

(0.3340) 

(0.4140) 

0:2 

0.4681 

0.3047 

0.1676 

0.1780 

0.0755 

0.3646 

(s.e.) 

(0.1879) 

(0.1825) 

(0.1792) 

(0.1717) 

(0.1044) 

(0.1659) 

ds 

1.1704 

0.7205 

0.4009 

0.4264 

0.1792 

0.8615 

(s.e.) 

(0.2866) 

(0.3811) 

(0.4053) 

(0.3843) 

(0.2443) 

(0.2929) 

A 

3.8994 

5.0600 

4.6565 

4.6252 

4.8877 

5.3466 

(s.e.) 

(0.5603) 

(0.8472) 

(0.6698) 

(0.6620) 

(0.7579) 

(1.0009) 

9 

— 

0.6932 

4.1980 

0.1484 

0.8212 

0.8844 

(s.e.) 

— 

(0.2318) 

(3.7768) 

(0.1351) 

(0.2739) 

(0.1250) 

log(^) 

-20.5904 

-18.6775 

-18.6864 

-18.7065 

-18.5941 

-18.4599 

AIG 

49.1809 

47.3549 

47.3728 

47.4130 

47.1882 

46.9198 

AICC 

51.0856 

49.2904 

49.3083 

49.3485 

49.1237 

48.8553 

BIG 

54.2133 

55.4095 

55.4274 

55.4676 

55.2428 

54.9744 

K-S (Ki) 

0.1034 

0.0952 

0.0865 

0.0871 

0.0837 

0.1011 

(p-value) 

(0.8240) 

(0.8906) 

(0.9448) 

(0.9418) 

(0.9576) 

(0.8435) 

K-S {Y 2 ) 

0.1001 

0.0900 

0.0957 

0.0961 

0.0920 

0.0843 

(p-value) 

(0.8527) 

(0.9255) 

(0.8870) 

(0.8844) 

(0.9132) 

(0.9551) 

K-S (max(Yi,y2) ) 

0.1431 

0.1405 

0.1329 

0.1331 

0.1344 

0.1518 

(p-value) 

(0.4344) 

(0.4587) 

(0.5301) 

(0.5289) 

(0.5161) 

(0.3613) 

LRT 

— 

3.8258 

3.8079 

3.7677 

3.9925 

4.2609 

(p-value) 

— 

(0.0504) 

(0.0510) 

(0.0522) 

(0.0457) 

(0.0389) 


For choosing the best value for k in BGEB distribution, we considered k = 2,..., 30, and 
obtained the corresponding log-likelihood values for all models. The results are given in Figure 
[2j It can be concluded that the log-likelihood values increase when k increases. But there is 




















negligible variation for large k. Therefore, the results for k = 30 are given in Table [2j 

Similarly, for choosing the best value for k in BGENB distribution, we consider k = 
1,..., 100. The results are given in Figure[2j It can be concluded that the largest log-likelihood 
value occurs in the case of = 2. Therefore, we present the results for k = 2 in Table [2l 

For each fitted model, the Akaike Information Criterion (AIC), the corrected Akaike in¬ 
formation criterion (AICC) and the Bayesian information criterion (BIC) are calculated. We 
also obtain the Kolmogorov-Smirnov (K-S) distances between the fitted distribution and the 
empirical distribution function and the corresponding p-values (in brackets) for Yi, Y 2 and 

max(li,l2)- 

Finally, we make use the likelihood ratio test (LRT) for testing the BGF against other 
models. The statistics and the corresponding p-values are given in Table [2l 
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